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A GENERALIZED CORRELATION OF VAPORIZATION TIMES OF DROPS IN FILM BOILING ON A FLAT PLATE 
by Kenneth J. Baumeister,* Thomas D. Hamill* and Glen J. Schoessow t 
Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 

ABSTRACT 

A theoretical dimensionless correlation for the vaporization times of discrete liquid masses in the 
Leidenfrost state is presented and verified experimentally with data in the literature. The correlation 
is presented as a single curve relating a dimensionless vaporization time to a dimensionless initial 
liquid volume. The correlation works well for the entire gamut of initial liquid volumes from spherical 
drops to large pancaked blobs for which experimental data exist. 


AUSZUG 

Eine theoretische dimensionslose Beziehung fur die Verdunstungszeiten diskreter flussiger Massen im 
Leidenfrostzustand, die im Versuch mit Angaben im Sehrifttum uberpruft wurde, wird angegeben. Die 
Beziehung zeigt sich als eine Einzelkurve des Verhaltnisses zwischen einer dimensionslosen Verdunstunggzeit 
und einem dimensionslosem flussigem Anfangsvo lumen. Die Beziehung ist anwendbar auf dem ganzen Bereich 
flussiger Anfangsvolumen, von spharischen Tropfen bis zu flachkugeligen Flecken, fur die experimentelle 
Ergebnisse vorhanden sind. 


AHHOTAUHH 

npnBO,n,HTCH h noATBepscflaeTCfl; 3 k c ne pHMe h t a jl hhmh aslhhhmh b JiHTepaType , TeopeTH^ecKan 
6e3pa3MepHaa KoppejiHi^HH BpeMeH ncnapeHUH ^HCKpeTHbix jkhakhx Macc b coctohhhh JIeK,neH$pocTa. 
Koppejiai^HH npe^cTaBJiHeTCH KaK eAUHan KpHBaa CBH3biBaiomaH 6e3pa3MepHoe BpeMH ncnapeHHH c 
6e3pa3MepHHM Hawaii bhbim ofi'beMOM sch^kocth. KoppejinuiHH cnpaBefljiHBa ajih Bceft raMMbi Ha^ajibHbix 
oS-BeMOB stHflKocTH ot c$epnwecKnx Kanejib .no 6ojibmHX npnnjiiocHyTbix Kanejib £Jih KOTopbix HMeioTCH 
3KcnepHMeHTajibHbie flaHHbie. 


INTRODUCTION 

If a quantity of liquid is placed on a suffi- 
ciently hot plate, the liquid will evaporate in 
the immediate vicinity of the plate at a rate 
sufficient to support the liquid. This phenomenon 
is referred to in the literature as either Leiden- 
frost boiling or film boiling. There are two 
broad categories of general interest that experi- 
mental work falls into ( Fig . 1 ) : 

(1) A continuous or infinite amount of liquid 
resting on the plate 

(2) Discrete or finite amounts (drops) of 
liquid. 

The discrete range (Figs. 1(a) to (e)) has 
many interesting problems associated with it that 
are absent from the continuous range. Figure 1 
shows a series of possible states that belong to 
the discrete and the continuous ranges. The addi- 
tional important variables in the discrete range, 
which are not associated with the continuous 
range, are the volume of the liquid mass and the 
drop shape; that is, the experimentalist working 
in the discrete range has one additional independ- 
ent variable, the volume of the liquid placed on 
the plate. For very small volumes (Fig. 1(a)) the 
shape of the drop is nearly spherical. With 
larger volumes (Fig. 1(b)) the drop tends to flat- 
ten out into a disk. The thickness of large 
masses of liquid, called extended drops, tends 
toward an asymptotic value ( Fig . 1(c)). For very 
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large volumes, bubble breakthrough starts occur- 
ring (Figs. 1(d) and (e)). 

Many of the problems associated with the dis- 
crete boiling range are examined in [l to 11].** In 
particular, Patel and Bell [7] show that bubble 
breakthrough (Figs. l(d) and (e)) is consistent 
with the prediction of instability theory for sub- 
merged surface film boiling. A numerical proce- 
dure for calculating the vaporization times of 
spherical drops (Fig. 1(a)) is presented in [8], 
while a theoretical analysis of heat transfer to 
large and extended drops (Figs. 1(b) and (c)) is 
presented in [5] and [9]. Recently, Baumeister, 
Hendricks and Hamill showed that metastable 
Leidenfrost film boiling (Figs. l(a) to (c)) can 
occur for plate temperatures as low as the satu- 
ration temperature of the liquid [ll] . 

This paper presents a theoretical correla- 
tion along with experimental confirmation that 
allows an a priori prediction of the heat-transfer 
area, heat- transfer coefficient, and vaporization 
time for any fluid, wall temperature, and liquid 
volume . 


METHOD OF ANALYSIS 

The vaporization time of a discrete liquid 
drop in Leidenfrost boiling on a flat plate can be 
found by a direct integration of the interface 
energy balance: 

= h T (V)A(V)AT (1) 

The functional form of the heat- transfer coef- 
ficient will be that presented in [5 and 9] . In 
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these references , a cylindrical model is used to 
predict the heat transfer. From [9] 


h c = 


0.68 


^ 3 A*gp z p v \ 1//4 

V AT^L e j 


( 2 ) 


where the properties are evaluated at the film 
temperature. 


V / V 

L p “ ■■■ 0 0 = geometry factor (3) 

* 2 l 2 


and 



The heat- transfer expression given by Equa- 
tion (2) is an analytical expression for the heat- 
transfer coefficient that would exist if no radia- 
tion from the plate to the drop existed; however, 
radiation is an important heat -transfer mechanism 
to consider, especially for high plate tempera- 
tures. Therefore, the total heat-transfer coeffi- 
cient to be used in the integration of equation 

(l) is given by 


hp - — (5) 

where f is the radiation correction factor. The 
radiation correction factor f is presented in 
the Appendix. For convenience and mathematical 
simplicity. Equation (l) is integrated by using a 
conductive heat- transfer coefficient rather than 
the total heat-transfer coefficient. After Equa- 
tion (l) is integrated, a mean value of the radia- 
tive correction factor f is used to give the 
correct vaporization time. 

Equation (l) still cannot be integrated since 
the heat-transfer area A and effective geometry 
factor L in the heat-transfer expression 
(Eq. (2))' 'are unknowns . Thus, the problem of 
determining these quantities is of paramount impor- 
tance in computing the vaporization time. A de- 
tailed method for computing the shape of the liquid 
drop by solving the Laplace capillary equation (see 
Fig. (2)) is presented in [4]: 



total vaporization time t c which assumes 
only a conductive and convective heat- 
transfer mode to the drop. 

(5) Apply a mean value of the radiation cor- 
rection factor f (presented in the Appen- 
dix) to the expression for the conductive 
vaporization time to give the actual vapor- 
ization time. 


DROP SHAPE 

A liquid drop not wetting its supporting sur- 
face, as in Leidenfrost film boiling, has the gen- 
eral shape depicted in Figure 2. The shape can be 
found by solution of the Laplace capillary equation 
(Eq. (6)). Equation (6) was nondimens ionalized and 
transformed into a coupled pair of ordinary non- 
linear differential equations, which relate the 
radius of curvature and the pressure difference at 
any point on the surface of the drop. The coupled 
equations were solved by a numerical integration on 
an IBM 7094 digital computer. The details of the 
transformation and numerical integration of the 
equation can be found in [4] along with an experi- 
mental verification of the numerical results. 

For a given drop volume, the maximum drop ra- 
dius and an average drop thickness l , defined by 
the equation 


l 



(7) 


were determined from the numerical solutions. 

The numerical results are shown as a dashed 
curve in Figure 3. The dimensionless groups shown 
in Figure 3 are defined as follows: 



From the universal curve, the average drop 
thickness and maximum radius can be determined pro- 
vided the surface tension and density of the liquid 
are known. For purposes of calculation, this curve 
is broken into three asymptotic ranges. 


The Laplace capillary equation relates the 
surface tension forces and the pressure difference 
across the surface of a drop. From the drop shape, 
the geometric parameters A and Lg can be deter- 
mined. 

In summary, the vaporization time for a liquid 
drop of volume V can be calculated by integration 
of the interface energy balance given by Equation 
(l). The procedural steps necessary to perform 
this integration are as follows: 

(1) Solve Equation (6) to determine the shape 
of the liquid drop for a given liquid volume 
V in terms of the basic fluid properties. 

(2) Determine the heat-transfer area and L e . 

(3) Substitute the value of L e into the heat- 
transfer equation (Eq. (2)) thereby giving an 
expression for the conductive heat-transfer 
coefficient in terms of the basic fluid 
properties and the drop volume V. 

(4) Substitute the functional form of he and 
A into Equation (l) and integrate for the 


Extended Drop Domain 

The extended drop domain is defined as the 
domain in which the thickness of the drop is in- 
dependent of drop volume. For extended liquid 
drops the universal curve (Fig. 3) approaches a 
slope of 3, which agrees with the physical observa- 
tions that the thickness of the drop approaches a 
constant asymptotic value. From Figure 3 the ex- 
tended drop region is defined for the domain of V* 
where 

V* > 155 (10) 

The average drop thickness in this domain is given 
by 
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This equation is independent of the volume of the 
drop and provides a method for computing the asymp- 
totic drop thickness. 

In addition, according to step 2 in the sec- 
tion METHOD OF ANALYSIS, the heat -transfer area for 
large cylindrically shaped extended drops can be 
defined as 

A = } (12) 

Therefore, when Equation (11) is substituted into 
Equation (12), the heat-transfer area of an extended 
liquid drop is 



Large Drop Domain 

The large drop domain is defined as the domain 
where the drop is disk- like in shape and where the 
thickness varies as a function of its volume, as 
compared with the extended drop whose thickness 
remains constant. The large drop region exists in 
the volume range where 

0.8 < V* < 155 (14) 


1* = -JU- = 0.83 (19) 

As can be seen in Figure 3, the universal 
curve approaches this asymptotic value of 0.83. 

From Figure 3, the small drop domain is taken to 
exist in that region where 

V* < 0.8 (20) 

Thus, in this region the drop thickness is given by 
l = 0.83 V 1 / 3 (21) 

The heat-transfer coefficient (Eq. (2)) de- 
rived in [9] was based on the assumption that a 
uniform gap existed beneath the drop; however, this 
is clearly not the case when a small drop is rest- 
ing on the surface. Consequently, the effective 
heat -transfer area is expected to be greater than 
the projected area of the sphere but less than the 
surface area of the lower half of the sphere; that 
is, 

*r 2 < A < 2itr 2 (22) 

Therefore, as an engineering approximation, 
the effective heat -transfer area of the sphere was 
taken as the average of the projected and surface 
areas; that is, 


as can be seen in Figure 3. The average drop 
thickness in this region is 


itr + 2nr 


= 1.5 nr* 


(23) 


l = 0 .b(^s\ V 1 / 6 (15; 

Vp z e/ 

In the same manner as in the extended mass 
region, the heat -transfer area in this region is 
also defined by Equation (12). Substituting Equa- 
tion (15) into Equation (12) gives for the heat- 
transfer area 


The effect of the preceding approximation was 
tested against the experimental data of Gottfried 
[2] and is reported in a later section entitled 
UNIVERSAL VAPORIZATION TIME PROFILE. 

For a sphere, 


r 



(24) 


A = 1.250^ V 5 / 6 


(16) 


This equation is applicable for drops intermediate 
in size, between small spherical drops and extended 
drops of constant thickness. 


Small Drop Domain 

When small quantities of liquid are in film 
boiling on a hot plate, their shape is nearly 
spherical. Thus, the average drop thickness de- 
fined by Equation (7) goes to a limiting value for 
small spherical drops: 


4 3 

, V 3_Y_ 4 

A 2 3 r 

itr 


But for a sphere 


V l/3 = (f) 


1/3 


( 1 ?) 


(18) 


Therefore, dividing Equation (17) by Equation (18) 
gives for small drops 


Substituting Equation (24) into Equation (23) gives 

(25) 

as the effective heat-transfer area of a sphere. 


A = 1.5(^f 3 V*/3 


HEAT-TRANSFER COEFFICIENTS 


Extended Drop Domain 


Substituting the value of Z for large drops 
(Eq. (ll)) into Equations (2) and (3) gives 


L - 64 \ — ^v — ) 


(26) 


for the domain defined by Equation (10). Surpris- 
ingly, the heat- transfer coefficient predicted in 
this case is independent of gravity. Physically, 
therefore, the gap thickness under an extended 
liquid drop is independent of the gravitational 
environment. In a larger gravitational field, 
however, the drop thickness will be thinner and 
the heat-transfer area larger, thereby giving rise 
to shorter vaporization times in large gravita- 
tional fields. 
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Note that the preceding observation does not 
apply to the near zero-gravity case, since the drop 
would not obtain an equilibrium position above the 
plate because of the reactive force of the gener- 
ated vapor. 


V kV/Vp/ff 3 / 


( 32 ) 


Large Drop Domain 


Substituting the value of l for large drops 
and Equation (15) into Equations (2) and (3) gives 

/kW/2pl/2p al/Zgl/sX 1 / 4 

hc = ) (27) 

for the domain defined by Equation (14). 


Small Drop Domain 

Substituting the value of l for small drops, 
Equation (2l), into Equations (2) and (3) gives 


h c - 


1. 


, 3 * 

k A gP z P, 
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1/4 


(28) 


for the domain defined by Equation (20). Note that 
the heat-transfer coefficient in the small drop 
domain is independent of surface tension because of 
the fact that all drops if small enough are spheres 
independent of the surface tension. The geometry 
factor Lg appearing in the heat- transfer coeffi- 
cient is also independent of surface tension for 
spherical drops. 


VAPORIZATION TIMES 


Smal 1 Drops 

The total vaporization time of a liquid drop 
can be found by a direct integration of Equation 
(1). Substituting Equations (25) and (28) into 
Equation (1), integrating, and solving for the 
vaporization times give 


t 


c 



(29) 


for the conductive vaporization time. 

Thus, according to the definition given to the 
radiation correction factor of Equation (Al) in the 
Appendix, the actual total vaporization time is 
given by 


t = 1.21 


\ i7 ‘ 

: 3 g**p v ar 3 / 


fyS/12 


(30) 


Rewriting the Equation (30) in terms of V* defined 
by Equation (8) gives 


t = 1.21 


iY'Vp^ 




iv 


(31) 


If a dimensionless total vaporization time t* 
is defined as 


the dimensionless total vaporization time for the 
small drop region is given by the equation 


t = 1.21 V 


*5/12 


(33) 


for the domain defined by Equation (20). 


Large Drops 

Substituting Equations (16) and (27) into 
Equation (l), evaluating the lower integration 
limits with Equations (14) and (33), applying 
the radiative correction factor, and introducing 
the nondimens ional variables defined by Equa- 
tions (8) and (32) give 

t* = 2.23 V* 1 / 3 - 0.97 (34) 

for the domain defined by Equation ( 14 ) . 


Extended Drops 

Substituting Equations (13) and (26) into 
Equation (l), evaluating the lower integration 
limits with Equations (10) and (34), applying the 
radiation correction factor, and introducing the 
dimensionless variables defined by Equations (8) 
and (32) give 

t* = 4.52 V* 1 / 4 - 5 (35) 

for the domain defined by Equation (10). 


Real Time Plots 


Figures 4 and 5 present a comparison of the 
theoretical and experimental total vaporization 
times for benzene in terms of real time and volume, 
as calculated from the equations presented in the 
previous sections. 

Theory and experiment are in good agreement 
for small drops (Fig. 4), for large and extended 
drops (Fig. 5), and even for drops with bubble 
breakthrough (10-ml drop). 


UNIVERSAL VAPORIZATION TIME PROFILE 

A plot of t* as a function of V* over the 
complete range of independent variable V* is 
shown in Figure 6. Some of the data in [1,2,4, 
and 6] have been plotted in terms of the dimension- 
less variables. The universal relation correlates 
quite well the available vaporization data over the 
entire gamut of liquid volumes. 

Note that no experimental data found in the 
literature fall in the small drop range defined by 
V* < 0.8; however, Gottfried’s experimental data 
came very close to falling in this range. 

Even though none of the initial liquid volumes 
reported in the literature fall into the small drop 
regime (V* <0.8) much of Gottfried's data is on 
the borderline of this limit (1.5 < V* < 5.5). 
Calculations show that for the smallest V* 
reported by Gottfried, the drop spends roughly 
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75 percent of its life in the small drop domain. 

The computed vaporization times for Gottfried’s 
data are, therefore, very dependent on the value of 
the heat- transfer coefficient and area computed for 
the small drop domain. The verification of the 
heat- transfer area approximation (Eq. (23)) is 
based on the agreement of calculated data with 
Gottfried's data, as shown in Figure 4. 

Two empirical correlations derived from dimen- 
sional analysis for the vaporization times of small 
drops are discussed in [2 and 8]. In one corre- 
lation Fick's law diffusion coefficient appears 
however, the diffusional dependency predicts an 
infinite evaporation time for experiments conducted 
in a saturated vapor atmosphere where there is no 
diffusion [ 12 , p. 376 and 4, p. 8]. This is con- 
trary to fact. The other empirical correlation 
predicts trends in agreement with the present cor- 
relation (Eq. (29)), although the functional de- 
pendency of the parameters is different. 


DISCUSSION OF RESULTS AND CQNC HIDING REMARKS 

The generalized correlation (Fig. 6) for pre- 
diction of vaporization times of discrete liquid 
masses presented in this paper is in good agreement 
with experimental data. Moreover, the validity of 
the correlation with respect to the volume of the 
discrete liquid drops placed on the hot surface 
applies over the whole spectrum of discrete states 
observed experimentally (Figs. 1(a) to (e)). 

Bubble breakthrough as shown in Figures l(d) 
and (e) seems to have a relatively minor effect on 
heat transfer. On the average, a blob of liquid in 
the extended drop region with bubble breakthrough 
appears to be equivalent to a flat disk. In all 
probability, bubble breakthrough does not alter the 
heat-transfer area. The presence of holes merely 
increases the perimeter of the bubbly liquid mass 
while the area remains constant. If the gap thick- 
ness beneath the drop remains fairly constant, the 
heat-transfer coefficient will not be altered very 
much. The net result appears to be that the total 
flux of heat input to the bubbly drop is nearly 
equal to that calculated by assuming no bubble 
breakthrough. This is not the case for a confined 
liquid (Fig. 1(f)) since bubble domes decrease the 
effective heat -transfer area. 

The vertical scatter of the data points shown 
in Figure 6 appears to result from the experimental 
uncertainties in the measured volume, vaporization 
time, and surface temperature. Most important, the 
surface temperature is reported as a constant value 
even though the actual plate temperature varies be- 
cause the drop cools the plate [6, p. 33]. There- 
fore, the experimental spread in the actual data 
(see Figs. 4 and 5) will be, of course, mirrored in 
the correlation shown in Figure 6. 

The temperature functionality in the universal 
vaporization time relations adequately correlated 
the data as seen by the real time plots (Figs. 4 
and 5); consequently, inverting the correlation 
(i.e., computing surface temperatures from observed 
evaporation times) should give fairly accurate 
results . 

The good agreement with theory tends to sub- 
stantiate the use of heat -transfer areas obtained 
by solving the Iaplace capillary equation for an 
isothermal drop. 


NOMENCLATURE 
A area, cm 2 

C specific heat of vapor at constant pres- 

sure, cal/(g)(°K) 

f radiation correction factor 

g acceleration of gravity, cm/sec 2 

g c conversion factor, 

l(g)(cm)/( dyne) (sec 2 ) 
h heat-transfer coefficient, 

cal/( sec ) ( cm 2 ) ( °K) 

h c calculated conductive heat-transfer coef- 

_ ficient, cal/(sec)(cm 2 )(°K) 

h c average value of h c during drop lifetime, 

cal/ ( sec ) ( cm 2 ) ( °K) 

h^ actual conductive heat- transfer coeffi- 

cient, cal/(sec)(cm 2 )(°K) 
radiative heat-transfer coefficient, 
cal/( sec ) ( cm 2 ) ( °K) 

hf total heat- transfer coefficient, 

cal/ ( sec ) ( cm 2 ) ( °K) 

k thermal conductivity of vapor, 

cal/(sec)(cm)(°K) 

Le equivalent geometry factor (see Eq. (3)), 

cm 

l average calculated drop thickness (see 

Eq. (7)), cm 

Z* dimensionless drop thickness (see Eq. (9)) 

P static pressure, dynes /cm 2 

AP pressure difference across surface, 

dynes /cm 2 
r radius, cm 

r max maximum radius of drop, cm 

r^rg radi^. of curvature of drop surface, cm 

T temperature, °K 

AT T^ - T s , o K 

t total vaporization time of drop, sec 

t* dimensionless total vaporization time of 

drop (see Eq. (32)) 

t c total vaporization time of drop considering 

conduction (see Eq. (29)) 

V drop volume, cu cm 

V Q initial volume of liquid placed on hot 

plate, cu cm 

V* dimensionless drop volume (see Eq. (8)) 

z distance from plate, cm 

6 gap thickness beneath drop, cm 

liquid emissivity 

A heat of vaporization, cal/g 

A* modified heat of vaporization (see 

Sq. (4)), cal/g 

4 absolute viscosity of vapor, poise 

p density, g/cm 3 

a surface tension, dynes /cm 

o Boltzmann constant, 4.876X10 -8 

k cal/(m 2 )(hr)(oK4) 

Subscripts : 

Z liquid 

s evaluated at saturation conditions 

v vapor 

w wall 
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APPENDIX - EFFECT OF RADIATIVE TRANSPORT ON VAPORIZATION TIMES 


The conductive total vaporization time t c 
given by Equation (29) has been calculated by 
assuming that all the heat is transported by con- 
duction and convection in the thin vapor film. Ra- 
diative transport has been neglected entirely. 

When radiation becomes appreciable at about 450° C , 
the theoretical vaporization times will be high. 

For vail temperatures up to about 550° C the 
following correction factor to t c is suggested: 
let t equal the vaporization time when radia- 
tion is present. Then, let f be defined by the 
equation 

t = ft c (Al) 


ation is not present. 


When radiation is present. 


- k 


8T 

5? 


/ he W 

z=6 


(A5) 


because of the nonlinear interaction of radiation 
and convection-conduction heat transport. Some 
quantity h^ exists that does satisfy the equa- 
tion 




(A6) 


where f is a factor less than 1 that accounts for 
radiation. 

The factor f is also equal to the mean value 
of the ratio of the heat- transfer coefficient ob- 
tained by neglecting radiation to the coefficient 
obtained by retaining the radiation term. By ex- 
tension of the analysis in [9, p. 13, Eq. (45)] the 
following equation can be shown to be valid in this 
temperature range: 


and if this quantity were known, then 

fcr = be + \ U7) 

However, since h^ is given by a simple formula, 
a functional relation of the form 

hr = f(V h c ) (as) 







(A2) 


where h r is the radiation heat- transfer coeffi- 
cient. For the flat disk geometry 


^ 'JV - 


(A3) 


The emissivity of the liquid in the tempera- 
ture range being considered can be taken con- 
veniently as equal to one if data are lacking. The 
transmissivity of liquid layers greater than a few 
millimeters in depth is effectively zero, and the 
reflectivity of liquids for wavelengths in the 
infrared is quite small [13, p. 371]. Hence, the 
emissivity is very nearly one. 

The plate emissivity does not enter the radi- 
ation equation because the ratio of the drop to 
plate area is effectively zero [2, p. 14]. The 
bottom surface of the drop "sees” the entire plate; 
hence, the view factor based on the undersurface 
of the drop is unity, which leads to_Equation (A3) . 

The heat -transfer coefficient should be 

evaluated at half the initial liquid volume under 
consideration because heat -transfer coefficient 
is a weak function of volume, and hence this 
procedure will give a reasonable estimate of the 
effect of radiation over the total life of the 
drop. 

Equation (A2) also points out that the heat- 
transfer coefficients are not simply additive. 

For small wall temperatures the equation linear- 
izes to 


hp = h c + 


0.75 hp 

\ 20 A } 


(A4) 


Note that 1^. is not the actual convection 
coefficient that exists physically. It is a 
quantity calculated with the assumption that radi- 


is more desirable. This is precisely the content 
of Equation (A4), which should not be misconstrued 
as a violation of the additivity of coefficients 
for parallel heat -transfer mechanisms. 
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Figure 1. - Film-boiling states of liquid masses. 



Figure 2. - Schematic of nonwetting liquid drop on flat surface. 
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Figure 4 - Droplet vaporization time against temperature difference for benzene. 
Initial drop volume, V 0 ■ 0. 0116 cm^ (data from Gottfried E2]). 
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Figure 5. - Liquid vaporization time against temperature difference for benzene 
(data from Patel and Bell [6 and 71). 
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Figure 6. - Universal total vaporization time curve. 



